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Abstract 

A Steiner star for a set P of n points in R d connects an arbitrary center point to all points of P, while 
a star connects a point p e P to the remaining n—1 points of P. All connections are realized by straight 
line segments. Fekete and Meijer showed that the minimum star is at most y/2 times longer than the 
minimum Steiner star for any finite point configuration in R d . The maximum ratio between them, over 
all finite point configurations in M. d , is called the star Steiner ratio in M. d . It is conjectured that this ratio 
is 4/-7T = 1.2732 ... in the plane and 4/3 = 1.3333 ... in three dimensions. Here we give upper bounds 
of 1.3631 in the plane, and 1.3833 in 3-space, thereby substantially improving recent upper bounds of 
1.3999, and V2 — 10~ 4 , respectively. Our results also imply improved bounds on the maximum ratios 
between the minimum star and the maximum matching in two and three dimensions. 

Our method exploits the connection with the classical problem of estimating the maximum sum of 
pairwise distances among n points on the unit sphere, first studied by Laszlo Fejes Toth. It is quite gen- 
eral and yields the first non-trivial estimates below y/2 on the star Steiner ratios in arbitrary dimensions. 
We show, however, that the star Steiner ratio in R d tends to v2> the upper bound given by Fekete and 
Meijer, as d goes to infinity. Our estimates on the star Steiner ratios are therefore much closer to the con- 
jectured values in higher dimensions! As it turns out, our estimates as well as the conjectured values of 
the Steiner ratios (in the limit, for n going to infinity) are related to the classical infinite Wallis product: 
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1 Introduction 

The study of minimum Steiner stars and minimum stars is motivated by applications in facility location and 
computational statistics EJHHUEl. The Weber point, also known as the Fermat-Toricelli point or Euclidean 
median, is the point of the space that minimizes the sum of distances to n given points in M, d . The problem 
of finding such a point can be asked in any metric space. It is known that even in the plane, the Weber point 
cannot be computed exactly, already forn > 5 0[1O]]. (For n = 3 and 4, resp., Torricelli and Fagnano gave 
algebraic solutions.) The Weber center can however be approximated with arbitrary precision (U|9l, mostly 
based on Weiszfeld's algorithm fP71 . The reader can find more information on this problem in ifTTTl . and in 
the recent paper of the the first two named authors lfl2Tl . 

The maximum ratio between the lengths of the minimum star and the minimum Steiner star, over all 
finite point configurations in R. d , is called the star Steiner ratio in R. d , denoted by pd- The same ratio for a 
specific value of n is denoted by pd{n). Obviously pd(n) < Pd, for each n. Fekete and Meijer [ 15] were the 
first to study the star Steiner ratio. They proved that pd < V2 holds for any dimension d. It is conjectured 
that p2 = 4/-7T = 1.2732 . . ., and p% = 4/3 = 1.3333 . . ., which are the limit ratios for a uniform mass 
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distribution on a circle, and surface of a sphere, respectively (in these cases the Weber center is the center 
of the circle, or sphere) ifTBTl . By exploiting these bounds, Fekete and Meijer also established bounds on the 
maximum ratio between the length of the minimum star and that of a maximum matching on a set of n 
points (n even) in two and three dimensions (d = 2, 3). 

In a recent paper, the upper bounds on the Steiner ratios for d = 2, 3 have been lowered to 1.3999 in the 
plane, and to \/2 — 1CT 4 in 3-space f!2l . The method of proof used there was not entirely satisfying, as it 
involved heavy use of linear programming. Besides that, the proofs were quite involved and the improve- 
ments were rather small, particularly in three-space. It was also shown in lfT2l that the bound 4/-7T holds in 
two special cases, corresponding to the lower bound construction; details at the end of Section [2] 

In this paper we get closer to the core of the problem, obtain substantially better upper bounds, and 
moreover replace the use of linear programming by precise and much shorter mathematical proofs. Here 
we prove that p2 < 1.3631, and ps < 1.3833. Based on these estimates, we can then further improve the 
estimates given on r]2 and 7/3, using the method developed by Fekete and Meijer [15]. Our improvements are 
summarized in TableQ] Here min S denotes a minimum star and its length (with a slight abuse of notation), 
SS* denotes a minimum Steiner star and its length, and max M denotes a maximum length matching and 
its length (for an even number of points). Finally, our method yields the first non-trivial estimates below 
\/2 on the star Steiner ratios in arbitrary dimensions. Among others, we show that the upper bound \[2 on 
the star Steiner ratio given Fekete and Meijer is in fact a very good approximation of this ratio for higher 
dimensions d; thus in this sense the problem in the plane is the most interesting one. 



Ratio 


Lower bound 


Old upper bound 


New upper bound 


P2 


(min S) 1 SS* 


| = 1.2732... 


1.3999 


1.3631 f 


P3 


(min S)/SS* 


J = 1.3333. . . 


y/2- 10" 4 = 1.4141... 


1.3833 f 


Pi 


(min S)/SS* 


^ = 1.3581... f 


y/2 = 1.4142... 


1.3923 f 


P5 


(min S)/SS* 


f = 1.3714... f 


y/2 = 1.4142... 


1.3973 f 


pioo : (min S)/SS* 


1.4124... f 


y/2 = 1.4142... 


1.4135 f 


rj2 : min S/ max M 


% = 1.3333. . . 


1.6165 


1.5739 f 


r/3 : min S/ max M 


| = 1.5 


1.9999 


1.9562 f 



Table 1: Lower and upper bounds on star Steiner ratios (p2, P3, P4, P5, P100X an d matching ratios (7/2, 773) for some 
small values of d. Those marked with f are new. 

If qo, qi, . . . , q n -i are n variable points on the unit (radius) sphere in W 1 , let G(d, n) denote the maxi- 
mum value of the function Yl,i<j i- e -> tne tne maximum value of the sum of pairwise distances among 
the points. It was shown by Fejes Tdth [14], and rediscovered in lfT2l . that G(2, n) has a nice expression in 
closed form: 

However, only the simpler inequality G(2, n) < ^n 2 will be needed here. The exact determination 
of G(d,n) for d > 3 is considered to be a difficult geometric discrepancy problem [3, pp. 298], however 
estimates of the form G(d, n) < Cdn 2 /2 are known O, where a is the "constant of uniform distribution" 
for the sphere (H [3] 01 : a equals the average inter-point distance for a uniform mass distribution on the 
surface of the unit sphere in R d . In particular for d = 3, Alexander U] [2]] has shown that 

^r? - lfW/2 < G(3, n) < \ 2 - ~. (2) 

Here C2 = 4/7T, and C3 = 4/3. The connection with this problem is explained in the next section. 
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2 Stars in the plane 



Fix an arbitrary coordinate system. Let P = {pq, . . . ,p n -i} be a set of n points in the Euclidean plane, and 
let pi = (xi, yi), for % = 0, . . . , n — 1. Let SS* be a minimal Steiner star for P, and assume that its center 
c = (x, y) is not an element of P. As noted in Ifl5l . the minimality of the Steiner star implies that the sum 
of the unit vectors rooted at c and oriented to the points vanishes, i.e., YH=o ]fpr = 0- For completeness, 
we include here the brief argument (omitted in [ 15 ]). The length of the star centered at an arbitrary point 
(x,y) is 



L(x,y) 

If (x,y) is the Weber center, we have 



n-1 

E 

i=0 



a/ \x - Xi) 2 + {y - yi) 2 . 



ox Oy 



0. 



The two equations give 

n-1 



E 



and 



n-1 

E 



y-Vi 



yj{x- Xi) 2 + {y - yif 



0. 



(3) 



i=Q V( x - x i) 2 + (y- yd 2 i= o 

Our setup is as follows. Refer to Figure [T] We may assume w.l.o.g. that the Weber center is the origin 

= (0, 0). We may also assume that the Weber center is not in P, since otherwise the ratio is 1. We can 
assume w.l.o.g. that the closest point in P to o is po = (1,0), hence SS* = (1 + 5)n, for some 5 > 0. Let C 
be the unit radius circle centered at o. We denote by v a vector v, and by 11^1 its length. Write rt = opi, for 

1 = 0, 1, . . . , n — 1, and let oqi = r?Sh be the corresponding unit vector; i.e., qi is the intersection between 

\°Pi\ 

J 7 } and the the unit circle C. Let = | r^|, hi = \poq\\, and a[ = \popl \ , for i = 0, 1, . . . , n — 1. We have 
SS* = J27=o a i- Finally* let aj G [0, 2tt) be the angle between the positive x-axis and r^. 





Figure 1: Left: estimating the length of a star centered atpo = (1, 0) € P. Right: estimating pairwise distances. 
Henceforth © can be rewritten in the more convenient form 



n-1 



n-1 



^^cosaj=0 and sin a,j = 0. 



(4) 



Since the orientation of the coordinate system was chosen arbitrarily, such formulas hold for any other 
orientation. Note that this is also equivalent with sum of unit vectors vanishing: Y17=o °Qi = 0- Moreover, 
such formulas hold for any dimension d > 2. 

Let Si be the star (and its length) centered at p^ for i = 0, 1, . . . , n — 1, and let min S = mino<j< n -i Si 
denote the minimum star. Using the local optimality condition ©, Fekete and Meijer |[T5l show that if one 
moves the center of SS* from the Weber center to a closest point of P, the sum of distances increases by a 
factor of at most y/2; this bound which is best possible for this method implies that for any d > 2 (see lfT5ll 
for details): 

p d {n) < y/2, thus p d < y/2. (5) 

From the opposite direction, by considering a uniform mass distribution on the surface of a unit sphere in 
W 1 , one has for any d > 2: 

Pd>c d - (6) 

Our new argument is a nutshell is as follows. If 5 is large, we consider the star centered at a point in P 
closest to the Weber center, as a good candidate for approximating the minimum star. If 5 is small, we use an 
averaging argument to upper bound the length of the minimum star. In the end we balance the two estimates 
obtained. Applying the averaging argument (for small 5) leads naturally to the problem of maximizing the 
sum of pairwise distances among n points on the surface of the unit sphere (or unit circle). 

Theorem 1 The star Steiner ratio in the plane is less than 1.3631. More precisely: 

4 2V2-± 

-<P2< < L363L 

Proof. Consider an n-element point set P, and the previous setup. It is enough to show that 

minS" ^ 2y/2 — — 



SS* - 1 + y/2- - ' 

By the triangle inequality (see also Fig. [Uleft) ), we have < 6j + ctj — 1, for i = 0, . . . , n — 1. Hence 

n— 1 n—1 n—1 

So = Yj a 'i ^ 12( bi + ai-l) = Y,h + Sn. 

i=0 i=0 i=0 

By Lemma 4 in lfl31. the local optimality condition (0]) implies YllZo k — n V^- It follows then that 
5*0 < (\/2 + 5)n. Hence the star Steiner ratio is at most 

So <V2+l_ 



SS* ~ 1 + 5 

Since the local optimality condition X^o 1 cos ai = ^ n °lds for any d, © also holds for any d. Observe 

that 

lim ^— = y/2, and lim ^— = 1, 

^0 1 + 5 <5^oo 1 + 5 

with the above expression being a decreasing function of 5 for 5 > 0. 

Clearly, the sum of the lengths of the n stars (centered at each of the n points) equals twice the sum of 
pairwise distances among the points. 

n-l 

i=0 i<j 
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By the triangle inequality (see also Fig. [Bright) ) 

\PiPj\ < \PiQi\ + \<HQj\ + \ljPj\ = (a* - 1) + \QiQj\ + («j ~ !)■ 
By summing up over all pairs i < j, we get 

n—1 n—1 

< 2 J2 + 2 ( n - !) X^ a * - = 2 Z) + 2<5 ( n - 

i=0 i<j i=0 i<j 

Using the upper bound estimate on G(2, n) in (Q]) we obtain 

2n 4 

VSj< — + 25(n -l)n<-n 2 + 25n 2 . 

' tan 7?- vr 

i=0 2 " 

The minimum of the n stars, min 5, clearly satisfies: 



inin5 < ^ =q1 - < [- + 25 } n. 
n \ 7r 

It follows that the star Steiner ratio is at most 



minS $ + 25 

< - . (8) 

SS* ~ 1 + 5 w 



This estimate holds for any dimension d: the points q-i lie on the surface of the unit radius sphere B in 
centered at o (rather than on the unit circle C). Observe that 

| + 28 4 1 + 25 
lim & = -, and lim 2 = 2, 

<5^0 1 + 5 7T 5^oo 1 + 5 

with the above expression being an increasing function of 5 for 5 > 0. Therefore, by combining this 
observation with the previous observation following ([7]), we get 

minS ^ . / V2 + 5 | + 2 A 
—— — < max mm — — , + — . 

55* ~ <5>o I 1 + 5 1 + 5 J 

The maximum value is given by substituting for 5 the solution of the equation v2 + 5 = ^ + 25 (that is, 
by balancing the two upper estimates on the Steiner ratio given by inequalities |7]) and © ). The solution 
5 = ^2 - | = 0.1409 . . . yields 

, x 2V2-- 2^2 

p 2 (n) < =-^4- = 1-3630 . . . , thus also p 2 < =-2-r = 1.3630 .... (9) 

1 + V2-- 1 + V2-^ 

□ 

By the result in |Q3), 55* < 4= maxM. By our Theorem[H min 5 < ^^J^ ■ 55*. Combining the 
two upper bounds yields the following upper bound on r]2'. 

Corollary 1 The minimum star to maximum matching ratio in the plane (r)i) is less than 1.5739. That is, 
for any point set 

minS 2V2-| 2 
max M \ + y/2-± ^3 
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The best known lower bound for this ratio, is 4/3, see lfl~5l . 

According to lfl2l Theorem 2], the star Steiner ratio for a set of n points in the plane that lie on a circle 
centered at the Weber center is at most 

— 4 4 

2n . _ < _ 

tan^ vr vr' 

The same bound holds for any finite point set in the plane where the angles from the Weber center to the n 
points are uniformly distributed (that is, a, = 2iir/n, for i = 0, 1, . . . , n — 1) [1.2. Theorem 3]. We think 
that this is always the case, and thereby venture a slightly stronger version of the conjecture proposed by 
Fekete and Meijer [ 15 ] (however, this is for specific values of n): 

Conjecture 1 The star Steiner ratio for n points in the plane is 

— 4 

tan Wi 7T 



3 Stars in the space 

In this section we give estimates on the star Steiner ratio in 3-space (Theorem and in higher dimensions 
(Theorem [3]). 

Theorem 2 The star Steiner ratio in M 3 is less than 1.3833. More precisely: 

| <P3 < ^(16-3^2) < 1.3833. 

Proof. The method of proof and the setup is the same as in the planar case, so we omit the details. Let B be 
the unit radius sphere centered at o, analogous to the unit circle C. Now all the points q± lie on the surface 
of B. Using the upper bound estimate on G(3, n) in Q we get the analogue of Equation ®: 

^8 1 + 2^ 

SS* ~ 1 + 5 y ' 

Taking also ((T)) into account, we have 

. [V2 + S l + 25\ 
03 < max mm — — , 4 =- . 

' 6 - s>o y 1 + 5 1 + 5 J 

By balancing the two upper estimates in © and (TTOb as in the planar case yields 5 = y/2 - | = 0.0808 

and _ 

2y/2 - I Qyj2 - 4 2 . r~. 

P3< =-^r = -^= = -=(16- 3V2 = 1.3832.... (11) 

H 1 + V2- | 3^/2-1 17 V ' 

□ 

By the result in 03), SS* < x/2-maxM. By our Theorem|2 min S < ^(16-3^2) -SS*. Combining 
the two yields the following upper bound on 773 : 

Corollary 2 The minimum star to maximum matching ratio in 3-space (773) is less than 1.9562. That is, for 
any point set 

miniS 2 , 1—^1- 4 . /— 

< —(16 - 3V2)V2 = — (8V2 - 3) < 1.9562. 



maxM ~ 17 v ' 17 
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The best known lower bound for this ratio, is 3/2, see lfl~5l . 

The same method we used in proving Theorem Q] and Theorem |2l together with various approximations 
yield the following estimates on the star Steiner ratio in M. d : 

Theorem 3 Let 1 < c<j < 2 be the "constant of uniform distribution" for the sphere in R , d > 2. The star 
Steiner ratio in M rf is bounded as follows: 

Cd< Pd< 2 ^ Cd ? w /j ere Km Cd = lim p d = \/2. 

1 + V2 — Crf d^oo d^oo 

The following closed formula approximations hold: 

2^2- V2e~ 5(2^1) 



\/2e 4 < 2d - 3 ) < < 



1 + V2 - \/2e s^ry 



Proof. Note that by the same argument used in the proofs of Theorem Q] and Theorem [2] (equations (O and 
CD} ), we have 

2^2 -c d 

Cd< Pd< ~ 7= • (12) 

1 + V2 - c d 

In order to establish the limits, we start by computing the "constant of uniform distribution" c<j. Recall 
that Cd equals the average distance from a point on the unit sphere in R d to all the other points on the same 
sphere, for a uniform mass distribution. It is easy to verify that Cd is given by the following integral formula: 

2 f" /2 sm d - 2 (2a) - sin a da 

c d = — ro — • (13) 

/J r/2 sin d - 2 (2a)da 

Some initial values are 

vr 4 64 48 

ci = l, C2 = -7 = 1.2732..., c 3 = - = 1.3333... , c 4 = — = 1.3581 ... , c 5 = — = 1.3714 ... . 
4 3 157T 35 

In order to establish a recurrence on Cd, define 

rir/2 



Jo 



sin* a ■ cos- 7 a da, i, j > 0. 

Jo 

Some initial values are 

ir In 
°oo = 2' fl01 = = an = 2' ° 02 = ° 20 = 4' 

Expanding sin 2a yields then 

Q = • 

Recall that integration by parts leads to the well-known recurrence relations for djj, for i, j > 1: 

r /2 . i i J sin*" 1 a • cosJ +1 a *f* i-1 r' 2 . ,_ 2 , J 

a,,- = / sin a • cos- 7 a da = H / sin a • cos-' a da 

Jo » + 3 i + 3 Jo 



sin i+1 a • cos- 7 1 a 



« + 3 



H / sm a • cos- 7 a da. 

t + 3 Jo 



1 



Plugging these in the formula for immediately gives a recurrence for c^. For any d > 1: 

2 ' WI ' jfe=T • a rf-M-2 4ri 2 / 1 \ 

-23" • 25^2 • a d-2,d-2 4d 1 \ 4 « V 

Recall at this point the infinite Wallis product from number theory! 16 1 : 



tt _ -fr / I/.'" \ 2 2 ! ! (i (i 



2 11 \4A; 2 - 1/ 1 3 3 5 5 7 79 

fe=l 

Let 

4A; 2 



- / Ah 2 \ °° 

^=n(i^rr)' and n 

fe=l V 7 fc=n4 



4& 2 - 1 

+i x 

denote the partial finite and respectively partial infinite Wallis products, so that W n Z n = tt/2, for every 
n > 1. Our recurrence for Cd yields that a is an increasing sequence satisfying also 

Cd+iCd+2 = cic 2 W d , d>l. (14) 
Since is bounded, it converges to some limit c. The value of c can be obtained by solving the equation 



c 2 = c l c 2^ = 2. 

We thus have lim^oo Cd = \[2. Since Cd < Pd < v2, we also have lim^oo Pd = y2. From Equation 
(TBI ), we also get that for d > 3 



7T 



cic 2 W d _ 2 < (? d < cic 2 W d -i, or v / cic 2 W ( i_2 < c d < y^i^Wd-i- (15) 

Observe that 



^ 4/c 2 - 1 ~ 2 ^ 1 2fc - 1 2k + l) 

=n+l fc=n+l v 7 



2(2n + 1) 



fc=n+l 

Standard inequalitie^] e 4;r / 5 < 1 + x < e x for x G [0, 1/3] now imply that for each n > 1 



fc=n+l v 7 



oo 1 



ECO 
k=n+l 



+ 1 4fc 2 -l = g2(2n + l) 



and 



I. ... I \ / 



g 5(2n + l) 

\ 1A-- - I / ~ 

k=n+l 

Since W n = {ir/2)/Z n , we have 



7T 1 7T _^_2 

_ . g 2(2n + l) < W < _ . g 5(2n + l) 

2 ~ 2 



and consequently (TT5l) gives 



2 v/vr I 2 v/vr l 

—= ■ 2-= ■ e 4 < M - 3 ) < Cd < —= • • e s(2d^TT 
V^r ^2 ~ y/n y/2 



'Here we have chosen 4:r/5 for simplicity of resulting expressions. 
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or equivalently 

^/2 e -i(2d-3) < Cd < \/2e~ 5 ( 2d -V . 



Taking into account (fT2|) and subsituting the above upper bound on a, we finally get the estimate (for any 

d > 4): 

r i 2\/2 - \/2e^s(2d3T) 

V2e 4 < M - 3 ) < p d < — I • 

1 + V2- ^e~^d=r> □ 



For the values of given by (U31 . we can extend the conjecture of Fekete and Meijer to all dimensions 

d> 2: 

Conjecture 2 The star Steiner ratio in M. d equals the "constant of uniform distribution " for the sphere in 
M. d : that is, pd = a for any d > 2. 
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